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ABSTRACT 

Recent millimeter- VLBI observations of Sagittarius A* (Sgr A*) have, for the first time, directly 
probed distances comparable to the horizon scale of a black hole. This provides unprecedented access 
to the environment immediately around the horizon of an accreting black hole. We leverage both 
existing spectral and polarization measurements and our present understanding of accretion theory 
to produce a suite of generic radiatively inefficient accretion flow (RIAF) models of Sgr A*, which we 
then fit to these recent millimeter- VLBI observations. We find that if the accretion flow onto Sgr A* is 
well described by a RIAF model, the orientation and magnitude of the black hole's spin is constrained 
to a two-dimensional surface in the spin, inclination, position angle parameter space. For each of 
these we find the likeliest values and their 1-cr & 2-er errors to be a = o +0 - 4+0 - 7 , 6 



and £ 
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2 go , when the resulting probability distribution is marginalized over the others. 



The most probable combination is a = 0+°- 2+0 - 4 , 9 = 90°_ 4 o°_ 5 o° and £ = -14_ 7 o±\\o, though 
the uncertainties on these are very strongly correlated, and high probability configurations exist for 
a variety of inclination angles above 30° and spins below 0.99. Nevertheless, this demonstrates the 
ability millimeter- VLBI observations, even with only a few stations, to significantly constrain the 
properties of Sgr A*. 

Subject headings: black hole physics — Galaxy: center — techniques: intcrfcrometric — submillimeter 
— accretion, accretion disks 



1. INTRODUCTION 

Understanding the structure and dynamics of black 
hole accretion flows has remained a central problem in 
astrophysics. Only in the past decade have sufficient 
numerical resources existed to perform self-consistent, 
three-dimensional magneto-hydrodynamic (MHD) simu- 
lations capable of resolving the magneto-rotational insta- 
bility, believed to be responsible for angular-momentum 
transport in black-hole accretion flows. However, a true 
ab initio computation is still well beyond reach, result- 
ing in the need for a variety of simplifying assumptions 
(e.g., the suitability of MHD, importance of electron-ion 
coupling, properties of accelerated electrons, etc.). As a 
result, the number of applicable models has rapidly pro- 
liferated, many of which are capable of describing the 
variety of phenomena observed. This is due in large part 
to the inability of current observations to resolve horizon 
scales. Unfortunately, the compact nature of black holes 
makes it very difficult to access the inner-edge of black 
hole accretion flows. 

Sagittarius A* (Sgr A*), the bright radio point source 
coincident with the center of the Milky Way, is presently 
the best studied known black-hole candidate. Observa- 
tions of orbiting OB-stars, the closest of which passes 
within 45 AU of Sgr A*, have produced a measured 
mass of 4.5 ± 0.4 x IP * 5 and an Earth-Sgr A* dis- 
tance 4 of 8.4 ± 0.4 kpc (jGhez et al.ll2009t iGillessen etlll 
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4 The mass and distance measurements are strongly correlated, 
with mass scaling roughly as M oc D l s . 



l2008h . Already, there is stro ng evidence for the existence 
of horizon in this source dBroderick fc Naravanl 120061 : 
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iBroderick. Loeb fc Naravanl 12009ft . However, in many 
ways, Sgr A* is very different than the supermassive 
black holes in AGN. Unlike its active brethren, Sgr A* 
is vastly under-luminous, with a luminosity many orders 
of magnitude smaller than the Eddington luminosity. As 
a result, it is widely expected that Sgr A*'s accretion 
flow is quite different from those in AGN, though per- 
haps more indicative of the roughly 90% of supermassive 
black holes that are not presently in an active phase. 

Even when strong-gravitational lensing is accounted 
for, the apparent angular size of Sgr A*'s horizon is only 
55 ± 2 /ias, a factor of two larger than the next largest 
black hole (M87) and orders of magnitude larger than 
any other known black hole (including all stellar-mass 
black holes). Nevertheless, despite its tiny angular size, 
recent millimeter VLBI experiments ha ve successfully re- 
solved this scale (jDoeleman et al.ll2008| ) . Since only three 
telescopes were involved with these observations, the re- 
sulting u—v coverage of the measured visibilities is very 
sparse. As a result only two simple models of Sgr A*'s 
im age, a Gaussian and an annulus, were fit to the data 
bv lDoeleman et al.l (|2008ft . 

However, we have a great deal of additional informa- 
tion about Sgr A*, including its spectral and polarization 
properties. We may also require physical consistency in 
any model (which would likely rule out an annulus, for 
example). Furthermore, we would like to evaluate the 
ability of, and optimize for this purpose, future millime- 
ter and sub-millimeter VLBI experiments to constrain 
fundamental properties of the accretion flow in Sgr A*. 
This paper demonstrates the fitting and parameter esti- 
mation procedure for a simple radiatively inefficient ac- 
cretion disk model (RIAF). In particular, we show that 
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even from very sparse baseline coverage it is possible 
to robustly extract interesting parameters of a generic 
RIAF. A study that considers how future high frequency 
VLBI obse rvations can exte nd this work will be presented 
elsewhere (jFish et alJl2008t ). 

In the analysis presented he re, we make full use o f 
all the observations described in lDoeleman et all l)2008l ). 
These include two days of VLBI observations at 230GHz 
that targeted Sgr A* using an array consisting of the 
James Clerk Maxwell Telescope (JCMT) on Mauna Kea, 
the Arizona Radio Observatory Submillimeter Telescope 
(SMTO) on Mt. Graham in Arizona, and one 10, m dish 
of the Coordinated Array for Research in Millimeter- 
wave Astronomy (CARMA) in California. Robust de- 
tections of Sgr A* and correlated flux density mea- 
suremnets were obtained on the CARMA-SMTO and 
JCMT-SMTO baselines. No detections were found on 
the CARMA- JCMT baseline resulting in upper limits 
on the correlated flux density. In addition, contempo- 
raneous total flux density measurements at the same fre- 
quency were obtained by using the full CARMA array 
operating as a stand-alone instrument. Full details of 
the obse rvations, calibrat i on, an d data processing can be 
found in lDoeleman et al.l (|2008l ). 

In scction[2]we present the simple accretion flow model 
that we employ. Sections [3] and |4] discuss the data fitting 
procedure, including the Bayesian method by which we 
do the parameter estimation, and the generic constraints 
placed by the current VLBI results. Section[5]details how 
we define our uncertainties and presents our parameter 
estimates. Finally, we conclude in section [SJ 

2. RIAF VISIBILITY MODELING 

2.1. Accretion Flow Modeling 

Sgr A* transitions from an inverted, optically thick 
spectrum to a optically thin spectrum near millimeter 
wavelengths. This implies that Sgr A* is only becom- 
ing optically thin at 1.3 mm. Due to relativistic effects 
this tran sition does not occur isot ropically for orbiting 
gas (e.g.. iBroderick fe Loebl[2006al ), becoming optically 
thin on the receding side at longer wavelengths than the 
approaching side of the accretion flows orbit. As a con- 
sequence, the opacity of the underlying accretion flow is 
crucial to imaging Sgr A*'s accretion flow. 

Despite being diminutive in comparison to the Edding- 
ton luminosity for a 4.5 x 10 6 M Q black hole, Sgr A* is 
still considerably bright, emitting a bolomctric luminos- 
ity of approximately 10 36 erg/s. Thus, it has been widely 
accepted that Sgr A* must be accretion powered, imply- 
ing a minimum accretion rate of at least lO~ lo M0/yr. 
It is presently unclear how this emission is produced. 
This is evidence by the variety of models that have been 
proposed to explain the emissi o n characteristics of Sgr 
A* fe e.. iNaravan et all Il998t iBlandford fc Begelmanl 



1991 iFalcke fc Markofl 120001 I Yuan et all 120021 " l200a 
Loeb fc Waxmanl 120071) . Models in which the emis- 
sion arises directly from the accreting material have 
been subsumed into the general class of RIAFs, de- 
fined by the generally weak coupling between the elec- 
trons, which radiate rapidly, and the ions, which effi- 
ciently convert gravi tational potential energy into heat 
(jNaravan et al.lll998| ). This coupling may be sufficiently 
weak to allow accretion flows substantially in excess of 
the 10 -10 Mq/jt required to explain the observed lumi- 
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Fig. 1. — High-frequency radio data and the range of fitt ed model 
spectra. The red circles are taken from lYuan et a l. ( 2004) (and ref- 
erences therein), the yellow squares are from lMarrone and 
the blue triangle is the single-dish flux during the VLBI observa- 
tions discussed here. Since all of the data points were not taken 
coincidcntally errorbars on the radio data shown by the red cir- 
cles are indicative of the variability, not the intrinsic measurement 
errors (in contrast to the yellow squares which were determined 
by coincident observations, and the blue triangle not used in the 
fit) . The gray regions show the envelope of the model spectra, with 
the dark gray region showing a < 0.9, and the light gray showing 
0.9 < a < 0.998. Finally, the black line shows the spectrum for the 
spin and inclination shown in Fig. f2] 

nosity with a canonical radiative efficiency. 

Nevertheless, th e detection of polarization from Sgr 
A* ab ove 100 GHz JAitken et al.ll2000t iBower et al1l200lL 
l2003t iMarrone et al.1 l2006f ). and subse quent measure- 
ment s of the Faraday rota tion measure (jMacauart et al.l 
I2006t IMarrone etld1l2007f) , has implied that the accre- 
tion rate near the black hole is significantly less than the 
Bondi rate, requiring the existence larg e-scale outflows 
(|A^l200q iQuataert fc Gruzinovl 12000) . Therefore, in 
the absence of an unambiguous theory, we adopt a sim- 
ple, self-similar model for the underlying accretion flow 
which includes substanti al mass-loss. 

For concreteness, as in IBroderick fcToebl (|2006aD , we 
follow [Yuan et all (|2003[ ). and employ a model in which 
the accretion flow has a Keplerian velocity distribution, 
a population of thermal electrons with density and tem- 
perature 
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Fig. 2. — Left: An example image of the radiatively inefficient accretion flow around Sgr A* at 230 GHz, centered on the black hole, 
viewed at 60° from the black hole spin axis, which is oriented towards North (§ = 0), for a moderately rotating (a = 0.5) black hole. The 
image appears as an crescent due primarily to the Doppler beaming and shifting associated with the relativistic orbital velocities. The 
color-scheme is normalized such that red is bright, blue corresponds to vanishing intensity and the total flux from the image is 2.4 Jy. 
Middle: Visibility magnitudes associated with the image shown on top. For reference, the positions of the observations in the u-v plane are 
also show by the white points (the circle, diamonds, & squares correspond to the single-dish, SMTO-CARMA, & JCMT-SMTO detections, 
and the triangle corresponds to the CARMA-JCMT upper limit). Right: Scatter broadened visibility magnitudes. Again the u-v positions 
of the VLBI observations are shown. 

fit the average radio, sub-millimeter and near-infrared 
spectrum of Sgr A*. Ho wever, our procedure is dif fer- 
ent from that employed in iBroderick fe Loebl (|2006al) in 
two respects. Firstly, we keep the radial index of the 
non-thermal electrons fixed for all models. Secondly, the 
fitting is performed systematically for a large number of 
positions in the inclination-spin parameter space, yield- 
ing a tabulated set of coefficient values. For every incli- 
nation and black hole spin presented here this was pos- 
sible with extraordinary accuracy (reduce \ 2 < f in all 
cases and < 0.2 for many), implying that this model 
is presently significantly under-constrained by the qui- 
escent spectrum alone 5 . Subsequently, we obtained the 
appropriate model parameters for arbitrary spins and in- 
clinations via a high-order polynomial interpolation. 

The range of spectra for the models we employ here 
is presented in Fig. [H together with the radio & sub- 
millimeter data that was used for the fitting procedure. 
The shown spectra deviate from the fitted spectra much 
more significantly than the scatter in the fitted spectra 
themselves. This is due almost entirely to interpolation 
process, and is most extreme at high-spins, where the 
coefficient values begin to change rapidly (and thus have 
larger interpolation errors). Nevertheless, these do not 
produce substantially different images at 1.3 mm from 
their tabulated counterparts, and it is important to note 
that the associated error in the underlying PJAF model 
serves only to favor very high spins (a > 0.99), since 
in these cases it results in a larger accretion-flow photo- 
sphere and, as we shall see, the small size of the photo- 
sphere of the high-spin models is the primary reason they 
are excluded. Hence, in reality, such models are more 
strongly disfavored than shown here. All of these models 
are also capable of producing the Faraday rotation mea- 
sures observed, and thus the polarimetric properties of 

5 Part of the reason for this is almost certainly the fact that many 
of the radio fluxes were measured during different observational 
epochs, and thus the flux uncertainties are indicative of the source 
variability, not the intrinsic measurement error. For two of the 
data points in Fig. [T] (the yellow squares), this is not the case, 
having been measured coincidentally, and thus these play a much 
more significant role in constraining the RIAF model parameters. 



and a toroidal magnetic field in approximate (f3 = fO) 
equipartition with the ions (which produce the majority 
of the pressure), i.e., 

B 2 ! m p c 2 r s 

In all of these, r s = 2GM/c 2 is the Schwarschild radius, 
p is the cylindrical radius and z is the vertical coordinate. 
Inside of the innermost-stable circular orbit (ISCO) we 
assume the gas is plunging upon ballistic trajectories. In 
all of these e xpressions the rad ial structure was taken 
directly from lYuan et all (|2003l ) and the vertical struc- 
ture was determined by assuming the disk height is com- 
parable to p. Given a choice for the coefficients and a 
radiative transfer model, images may then produced us- 
ing the fully-relativistic ray-tracing and rad i ative t rans- 
fer schemes descri b ed in IBroderick fc Loebl (|2006aD and 
IBroderick fc Loebl (|2006bD . 

The primary emission mechanism is synchrotron, aris- 
ing from both the thermal and non-thermal electrons. 
We model the emission from t he thermal e l ectron s us- 
ing the emissivity described in lYuan et all (|2003h . ap- 
propriatel y altered to account for re lativistic effects 
(see, e.g.. IBroderick fc Blandfo7dll2004l ). since we per- 
form polarized radiative transfer via the entire com- 
plement of Stokes parameters, we employ the polar- 
ization fraction for the r mal s ynchrotron as derived in 
iPetrosian fc McTiernanl (|l983l ) . In doing so we have im- 
plicitly assumed that the emission due to thermal elec- 
trons is isotropic, which while generally not the case is 
unlikely to change our result s significantly. F o r the non- 
thcrmal electrons wc follow Uones fc O'Del (|l977l ) for 
a power-law electron distribution, cutting the electron 
distribution off below a Lorentz factor of I0 2 and corre- 
sponding to a spectral index of ayisk = 1.25, both roughl y 
in agreement with the assumptions in lYuan et alj (|2003h . 
For both the thermal and non-thermal electrons the ab- 
sorption coefficients are be determined directly via Kir- 
choff's law. 

As in IBroderick & Loeb (2006a), to correct for the 
fact that lYuan et all (|2003h was a Newtonian study, the 
three coefficients (n° th , T® and n° nth ) were adjusted to 
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Fig. 3. — Effective reduced \ 2 f° r the best fit models in £ and 
Vb as a function of inclination and spin. For reference, white lines 
on the color-bar show the minimum \ 2 (solid), minimum \ 2 + 1 
(dashed) and minimum \ 2 + 2 (dot-dashed) . We were able to find 
a good fit for nearly all inclinations above 40° . Models with lower 
inclinations are too large at all £ to produce the observed flux on 
the JCMT-SMTO baseline. Conversely, models with large inclina- 
tions and high spins were over-predict the fluxes on the SMTO- 
CARMA and CARMA-JCMT baselines. 
Sgr A*. 

At the time t hat S gr A* was being monitored by 
iDoeleman et al.l (I2008T ) it exhibited an anomalously low 
1.3 mm flux of 2.4 ± 0.25 Jy, almost 40% below its av- 
erage value (shown by the blue triangle in Fig. [TJ). To 
appropriately account for this we chose a "minimal" pre- 
scription for changing the model, reducing all densities 
by a fixed factor (decreasing the magnetic field such that 
(3 is fixed) until the observed flux was reproduced. In a 
sense, this models the low flux as a low-mass accretion 
period. These lower-density models do have noticeably 
smaller, though qualitatively similar, images compared 
to those used in the spectral fits. 

An example image is shown in the left-panel of Fig. 
[2] for a moderate spin and viewing inclination. As seen 
in previous efforts to image relativistic accretion flows, 
the flow appears as a crescent assoc iated with the ap- 
proach ing side of the accretion disk (jBroderick fc Loebl 
l2T)06allbh . The non-negligible optical depth obscures the 
black hole "silhouette" on this side. The receding side is 
all but invisible as a consequence of the Doppler shift and 
Doppler beaming associated with the relativistic orbital 
motion in the inner-most portions of the accretion flow. 
Images for 9 > 90° are related by a reflection (across 
a line perpendicular to the projected spin axis) to im- 
ages with an inclination 180° — 9. However, due to the 
approximate reflection symmetry of the RIAF model im- 
ages, the constraints arising from spins pointed counter 
to and along the line of sight are nearly identical. As a 
result, for clarity we restrict ourselves to the latter range. 

Using this model we created a library of ideal- 
resolution images, each with a flux of 2.4 Jy, having dif- 
ferent spins and viewing inclinations spanning all possi- 
ble values. 

2.2. Interstellar Electron Scattering 

The existence of an interstellar scattering screen be- 
tween Earth and the Galactic center has been well known 



for some time now. This has been carefully character- 
ized empirically b y a number of auth ors; we use the re- 
cent model from iBower et all (|2006t) . In this the ob- 
served flux distribution is obtained by convolving the 
ideal flux with an anisotropic Gaussian scattering ker- 
nel. The anisotropic Gaussian is defined by the scattering 
widths along major and minor axes (both cx A 2 ) and the 
position a ngle of the m i nor ax is (which is independent of 
A). From IBower et aD l|2006h . the associated full- width, 
half-max for the major and minor axes are 

FWHM maj = 1.309 ^Y^) mas ( 5 ) 

FWHM min = 0.64 ( —\ mas (6) 
\1 cm J 

with the major axis oriented 78° East of North. In prac- 
tice, the broadening was done in the u-v plane, where it 
reduces to a multiplicative factor. 

2.3. Visibility Modeling 
Visibilities arc then defined in the normal way: 

V(u, v) = J J dx dy e -** i ( xu+uv)/x I(x, y) (7) 

where I(x, y) are the intensities on the image plane, with 
x aligned East- West and y North-South. Because chang- 
ing the position angle, £ (vanishing at North and increas- 
ing towards the East), corresponds simply to a coordi- 
nate rotation, the images were originally computed only 
for a single position angle, namely £ = 0. The visibili- 
ties associated with the ideal-resolution images were then 
computed, and subsequently rotated to the desired £. In 
practice, the visibilities were calculated via a Fast Four- 
rier Transform, which was padded sufficiently to resolve 
the shortest baseline (SMTO-CARMA). We also vary the 
total flux, or cquivalently Vb = V(0,0). For sufficiently 
different values a new image must be produced. How- 
ever, for small variations around our canonical value of 
2.4 Jy (within 0.5 Jy) simply rcnormalizing the visibilities 
produces a good approximation. Finally, the interstcllar- 
scatter broadening was effected by multiplying the visi- 
bilities in the u-v plane by the Fourier-transformed scat- 
tering kernel. 

Since Sgr A* was detected on only two of the three 
VLBI stations, there was insufficient observational data 
to determine the individual baseline visibility phases. 
Thus, henceforth, by visibility (and V) we shall, more 
properly, refer to the visibility magnitudes. These are 
necessarily a function of position in the u-v plane, black 
hole spin (a), spin inclination (9), position angle of the 
projected spin vector (£) and flux normalization (Vb). 
That is, for a particular realization of the RIAF model, 
we have computed V(u, v; a, 9, £, Vq). 

The ideal and broadened visibilities arc shown in the 
center and right-panels of Fig. [2 respectively. In both of 
these the positions of the observed visibilities are shown 
by the white points. As expected, associated with the 
narrow-axis of the image crescent is a correspondingly 
broad feature in the visibilities. On the maximum scale 
of interest (the JCMT-SMTO baseline) the power is not 
substantially reduced by the interstellar electron scatter- 
ing. At the same time, along the long-axis of the image 
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Fig. 4. — Probability densities of a given inclination and black hole spin for our canonical Sgr A* model. The top left, top right and 
bottom left panels show the contribution to p(a, 0) from the SMTO-CARMA, JCMT-SMTO and CARMA-JCMT baselines, after being 
marginalized over £ and at the most-likely Vo. The combined implied p(a,8), marginalized over £ and Vo is shown in the bottom right. 
(Islands of high probability in the right panels are an artifact of under sampling the probability peak in inclination.) In each, p(a,8) 
is normalized such that the average is unity, providing a clear sense of the significance of the variations in probability. Note that the 
color-scale is different in each panel. While the JCMT-SMTO measurement is clearly the most constraining, the other baselines arc critical 
to eliminating high-inclination, high-spin solutions. 



crescent the visibilities drop off rapidly. Thus, we may 
expect that, at the very least, the VLBI observations will 
constrain £. This turns out to be correct. 

3. BAYESIAN DATA ANALYSIS 

We follow a Bayesian scheme to compute the prob- 
ability that given the measured visibilities, a given set 
of model parameters (a,f?,£) are correct. This is, of 
course, predicated upon the assumption that our simplis- 
tic RIAF model is appropriate for Sgr A*. Indeed, this is 
the primary uncertainty in our reported constraints upon 
a, 9 and £. While we can address this model dependence 
somewhat by comparing the results from different black 
hole masses, at this point we can only hope that our 
results are characteristic of generic RIAF models. This 
is not unreasonable given that the primary physics re- 
sponsible for the structure of our images, the Keplerian 
velocity profile, is a common theme among RIAF models. 
Nevertheless, it remains to be proven. 

Assuming the observational errors are Gaussian, the 
probability of measuring a visibility V. given a particular 



set of model parameters (a, 9, £, Vo) is 
Pi(Vi\aMV ) = 



■ exp ■ 



2AV/ 



dVi 



(8) 

This is appropriate for detections, i.e., along the SMTO- 
CARMA and JCMT-SMTO baselines. However, for the 
non-detection associated with the CARMA-JCMT base- 
line this is not the case. As discused in Appendix [AJ 
the probability of a non-detection given a measurement 
threshold of Vi and intrinsic uncertainty of AVi and ex- 
pected value V(ui, Uj; a, 9, £, Vo) is 



Pi{<Vi\a,e,£,V ) 



Vi - V(ui,Vi;a,9,£,V ) 
V2AV 



(9) 



Therefore, the probability of observing the measured set 
of independent visibilities given a particular RIAF model 
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Fig. 5. — Top: Probability densities of a given inclination and black hole spin assuming that Ms gr A* = 2.25, 4.5 and 9.0 X 10 Mq 
(note that the best current estimate is 4.5 ± 0.4 X 10 6 Mq). Bottom: Probability densities of a given inclination and position angle for the 
same masses. Note that the orientation of the accretion disk, and hence the black hole spin, is strongly constrained for all cases to lay in 
a taco-shell-shaped 2-dimensional surface within the inclination-spin-position angle parameter space. Regardless of mass, lower spins and 
moderate inclinations are preferred by the recent VLBI observations. A position angle of 0° corresponds to the projected spin vector being 
oriented North-South, and increases towards East. 



IS 



P{{V i }\a,e,i,V Q )= WIM.&Vo) 



i=SC, JS 

x 



II Pj{<Vj\a,0,e,V ),(10) 

j=CJ 

where the first product is over the detections on the 
SMTO-CARMA {SO) and JCMT-SMTO (JS) base- 
lines, and the second is over the non-detections on the 
CARMA-JCMT (GJ). 

In order to assess the quality of the modeling at each 
spin and inclination, we appeal to the Bayesian Informa- 
tion Criterion (BIC). Specifically, 



BIC = -2 In 



PM}\a,0,S,V o ) 



k lni 



(11) 



where k = 4 is the number of parameters and n = 20 is 
the number of data points. In the absence of the non- 
detection, this reduces to the normal definition of % 2 , 
up to an additive constant, and thus comparing mod- 
els with similar numbers of parameters reduces to the 
normal \ 2 minimization. Since this latter statistic is 
commonly used, and its properties generally well known, 
solely for the purpose of assessing the quality of the fits 
wc define an effective x 2 : 



-2 In 



P({Vi}\a,e,^V ) 



In (l[2nAV? (12) 



BIC + const. 



Note that since only one of the visibility measurements is 
an upper limit, we may expect that this definition of x 2 
will behave very similarly to the standard definition in 



our case. The reduced effective x 2 is shown as a function 
of spin and inclination in Fig. [3] for the best fit Vo and 
£. For nearly all inclinations above 35° a good fit can be 
found (indeed, \ 2 ' jv ~ 0.4!). The exception is at very 
high spins and high inclinations (where the disk is edge 



While eq. (fT0| gives the probability that the observed 
visibilities come from a given model, what we would like 
to know is somewhat different: the probability of a set of 
model parameters given the observed visibilities. That is, 
we would like the probability density p(a,@,£,Vo\{Vi}). 
With an appropriate choice of priors on a, 9, £ and Vo, we 
may construct the desired probabilities via Bayes' theo- 
rem. As such, we now turn to the problem of choosing 
these priors. 

A natural choice for the prior upon 6 and £ comes from 
the assumption that Sgr A*'s spin orientation probabil- 
ity is isotropic, i.e., we have no other information regard- 
ing its direction. This results in p(8, £) = sin 9. In the 
absence of a complete theoretical understanding of the 
spin evolution of supermassive black holes, we choose the 
prior on a to be uniform, i.e. p(a) = 1. Finally, we set 
the prior on Vb to be uniform as well. Since the allowed 
range of variation in Vo is small, and the prior probabil- 
ity is expected to vary smoothly, this is not a significant 

6 Note that in this procedure we have ignored the statistics as- 
sociated with fitting the spectrum, choosing to break up these two 
components of the fitting process to stress the implications of the 
new millimeter VLBI observations. 
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oversight. Therefore, Bayes' theorem gives 
p(a,6,Z,V \{V t }) 

P({V t }\a, 6, g, V o )p(a)p(0, Qp(V ) 
J dad9d^dV P({Vi}\a, 9, £, V )p(a)p(6, £)p(V ) 

= P({Vi}\a,O,Z,Vo)8m0 

J dad9d^dV P({V}\a,9,^V )sm9 

This is necessarily a probability density in a four- 
dimensional parameter space, and thus is quite difficult 
to visualize directly. Furthermore, some of these param- 
eters are physically more interesting than others. There- 
fore, we construct a variety of marginalized probabilities 
from this for presentation and analysis. The most general 
is simply marginalized over Vq: 

PM,0 = J W p{a,6,t,Vo\{V$). (14) 

This probability distribution is shown explicitly in Fig. 
[6l We construct a pair of two-dimensional marginalized 
probability densities as well: 

p(a,0) = J dZW a p(a,e,t,V \{Vi}) (15) 

and 

p(0,O = J dadV 0P (a,9,t,V \{V t }) (16) 

These are plotted in Fig. 's HI and 03 Alternatively, we 
could choose the most likely values of either Vq or £. For 
some of the panels in Fig. HI we choose a hybrid prob- 
ability: marginalized over £ but the most likely in Vq. 
Finally, for the purpose of identifying the probability dis- 
tribution of each parameter separately, we also construct 
the marginalized one-dimensional probability densities: 

p(a) = J d0d£dFoKM,£,^olW}) (17) 

P (9) = J dadtdV p{a,9,i,Vo\{Vi}) (18) 

p(0= J dad9dV oP (a 1 9,^Vo\{V l }). (19) 

These are shown in Fig [7] 

4. THE NATURE OF THE VLBI CONSTRAINTS 

The relative importance of the three different baselines 
is shown in Fig. HI To make a direct comparison of the 
contributions from different baselines possible, the top- 
left, top-right and bottom-left panels of Fig. HI show the 
probability after setting Vq to the most likely value and 
marginalizing over £ (justified by the fact that the prob- 
ability distributions in £ are quite similar, while those in 
V were not). Both the SMTO-CARMA detections and 
the CARMA-JCMT upper limit exclude the high-spin, 
high-inclination portion of the parameter space. This is 
primarily because high-spin RIAF models result in im- 
ages that are sufficiently compact to substantially over- 
predict the flux observed on both of these baselines. Due 



to its considerably longer baseline length, the CARMA- 
JCMT non-detection is more constraining than the mul- 
tiple SMTO-CARMA detections, despite only represent- 
ing an upper-limit. 

As anticipated, most constraining are the visibilities 
measured on the JCMT-SMTO baseline. This baseline 
convincingly excludes low inclinations (face-on disks). 
This may be easily understood, qualitatively, in terms of 
the RIAF images themselves. As inclination increases, 
the RIAF image becomes more asymmetric, and increas- 
ingly dominated by a thin crescent (e.g., left-panel of 
Fig. |2|). It is generally the direction across the minor- 
axis of the crescent that has the shortest scale inten- 
sity variations, and consequently determines the long- 
baseline visibilities. Below a critical inclination, the cres- 
cent grows sufficiently fat that it is resolved out by the 
JCMT-SMTO baseline. 

At large inclinations, the crescent becomes too thin, 
requiring the JCMT-SMTO baseline to be oriented 
obliquely relative to the spin axis of Sgr A*. This re- 
stricts the available values of £ to an increasingly limited 
range, as may be seen explicitly in the bottom-center 
panel of Fig. Thus, while it is nearly always possible 
to obtain a satisfactory fit, the reduced range in position 
angle makes such a configuration unlikely. Hence we find 
that the long-baseline visibilities strongly constrain the 
possible values of inclination and spin to a narrow band 
near 9 ~ 50° and generally favoring moderate a. We note 
that the appearance of islands of high probability in the 
top- right panel of Fig. HI is almost certainly due to under 
sampling in 9 near this critical region, and not indicative 
of a bifurcation in the allowed parameter space. 

This behavior is clearly visible in the combined 
p(a,6), which closely resembles that from the JCMT- 
SMTO baseline alone. In this case, we have marginal- 
ized over Vq as well as £ (i.e., p(a,6) as defined by 
cq. [T5|) . though the difference had we chosen the most 
likely value of Vq is less than 10% everywhere. The 
shorter baselines have effectively removed a small high- 
spin, high-inclination island that persisted in the JCMT- 
SMTO probability distribution alone. Additionally, they 
have further restricted the spin to low-to-moderate lev- 
els. The combined result, however, is to restrict the spin 
and inclination to a narrow strip. 

In a similar fashion, p{9, £) appears to limit RIAF mod- 
els to a narrow band of orientations. Due to the approx- 
imate up-down symmetry of the image (parallel to the 
projected spin axis), the corresponding band has a char- 
acteristic "U" shape, corresponding to something akin 
to a taco shell in the a,9,£ parameter space. (Note that 
we do not plot £ over the entire 360° range due to the 
symmetry of the visibilities under reflection.) 

In order to ascertain how the present uncertainty in 
the mass of Sgr A* effects our conclusions, we repeated 
the analysis for black hole masses of 4.1 x 10 6 Mq and 
4.9 x 10 6 M Q . Since the mass uncertainty is strongly 
correlated with the distance to Sgr A* (MD~ 18 is 
very well determined by stellar orbits. iGhez et aDl2009l : 
IGillessen et all l200l . we altered the distance to the 
Galactic center accordingly. As a consequence, the 
roughly 10% change in the mass results in a 4% change 
in the angular scale of the RIAF images. A correspond- 
ingly small change was seen in the resulting probability 
distributions, implying that, apart from the the RIAF 
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Fig. 6. — p(a, 8,£) is shown for slices of constant £. The slice with the maximum over-all probability is shown in the center panel. The 
white solid and dashed lines show the 1-cr and 2-cr contours, defined explicitly in the text, respectively. Note that while the probability is 
quite peaked, likely values of the parameters exist for a large number of spins, inclinations and position angles. 

modeling, the paucity of millimeter- VLBI observations 
is the dominant obstacle to constraining the black hole 
spin properties. 

Altering the angular scale of the images also provides 
a proxy, albeit a poor one, for considering different Sgr 
A* models, corresponding to different disk scale lengths. 
p(a, 9) and p(9, £) for a black holes of mass 2.25 x I 6 M Q 
and 9.0 x 10 6 M Q are compared to those determined us- 
ing the estimated value of 4.5 x f 6 Af© in Fig. [5] These 
extreme mass changes correspond to a 30% change in 
the angular scales of the images, and can make a sub- 
stantial difference to the resulting probability distribu- 
tions. Nevertheless, high inclinations and low spins are 
still generally preferred, though less so for large images. 

Finally, because Sgr A* is an inherently dynamical en- 
vironment, exhibiting substantial variations in the mil- 
limeter flux on 30 min time scales, we may not be jus- 
tified in assuming that the image of Sgr A* was sta- 
tionary during the entire time that observations were 
made. Indeed, searching for variations in the VLBI clo- 
sure quantities has been suggested as a way to directly 
probe the existence of hot spots in Sgr A*'s accretion 
flow (jDoeleman et al.l 12008). Therefore, to check this, 
we repeated this analysis for the two days over which 



the observations were performed separately, finding no 
significant difference. This implies that either Sgr A* was 
quiescent during this time or that the individual 3.5 hr 
observation windows were sufficiently long to average out 
this activity. The anomalously low 1.3 mm flux suggests 
the former interpretation is correct. 

5. PARAMETER ESTIMATION 

It is evident from the previous section that the allowed 
parameter space is highly non-Gaussian. As such, we 
must take special care in how we extract values and their 
attendent uncertainties for the fitting parmctcrs. In all 
cases, these values will be highly correlated, and the sys- 
tematic uncertainties due to the choice of a particular 
RIAF model will dominate the errors (and thus will not 
be reviewed again in this section). 

In order to determine l-cr and 2-cr error surfaces, we 
first define a cumulative probability 



P(>p) 



p(x) dx 



(20) 



p(X)>p 



where x are the parameters of the total probability dis- 
tribution. That is, P(> p) is simply the probability 
associated with the region of the parameter space that 
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Fig. 7. — From left to right: p(a), p(0) and all marginalized 

highly non-Gaussian. The dark and light shadded regions denote the 
to keep in mind that parameter estimates from these distributions ai 

has probability density above p. This is necessarily a 
monotonic function of p, which may then be inverted 
to find p as a function of P{> p). We then define the 
1-er contour to be that associated with the p for which 
P(> p) = 0.683. Similarly, we define the 2-cr contour by 
P(> p) — 0.954. While these satisfy the normal defini- 
tion of l-<7 and 2-cr errors, it is important to remember 
that in our case the errors are strongly non-Gaussian. 

Fig. [6] shows the probability density in the full three- 
dimensional parameter space as a sequence of constant- 
£ slices. The most likely parameter combination is 

a = +0 - 2+0 - 4 , 9 = 90°_4o°-50° and £ = -U t 7 7 l±n<,. 
However, as remarked in the previous section, there are 
acceptable solutions for a wide range of a, 9 and £. The 
primary effect of the VLBI measurements is to restrict 
these parameters to a band around a two-dimensional 
surface, resembling a taco shell with its sides around 
£ ~ —20° and £ ~ 63°, the highest probability densities 
clustered on the former. Generally, we see a preference 
for low spins and can rule out inclinations below 30°. 

More diagnostic for single parameters are the fully 
marginalized probability distributions shown in Fig. [71 
It must be remembered, however, that these parameter 
estimates are strongly correlated, and thus these esti- 
mates should be used with caution. Again we use the 
cumulative probability, P(> p) to define the 1-er and 2-cr 
surfaces, shown in the panels of Fig. [Jj as the dark and 
light shaded regions, respectively. 

From the left panel of Fig. [7J the most likely spin 
is a = o +0 - 4+0 - 7 . While we can rule out very high spins 
(a > 0.99) at the 2-cr level for our particular RIAF model, 
the spin is otherwise weakly constrained. Moreover, it 
appears to have a small high-spin island, at the 1-er level, 
around a = 0.9. The presence of this island is limited by 
the CARMA-JCMT non-detection, and therefore can be 
directly probed via future observations. 

In contrast, the inclination is fairly robustly limited. 
The most likely inclination is 9 = 50°j^o°iio° , where 
the lopsided errors are due to the lopsided nature of the 
probability distribution. It is very clear in the center 
panel of Fig. [JJ that face-on geometries (9 < 30°) are 
convincingly ruled out. However, there is a tail extending 



over all other parameters. In all cases the probability distribution is 
1-cr and 2-<r regions, as defined in the text, respectively. It important 
e strongly correlated. 

to higher inclinations 7 . 

The distribution of position angles is more complicated 
than either spin or inclination. In the right panel of 
Fig. [Jj we see the taco-shell geometry most clearly in 
terms of the bimodal distribution of likely £. The most 

I Q 1 O i 1 nwO 

probable position angle is £ = — 20_ 16 o_ 2 9° • However, 
a second solution does exist at £ = 63 _^ 7 o_ 112 o , though 
containing only 38% of the probability (under the Ter 
peak) of the more likely solution. 

Our limits upon the inclinati on are in quit e good agree- 
ment with previous efforts. iHuang et al.l (|2007f ) fit a 
qualitatively different RIAF model, primarily neglecting 
orbital motion, at 3.5 mm and 7 mm. While these favor 
high inclinations, the fact that the images are dominated 
by interstellar scattering prevents a conclusive measure- 
ment. Similarly, efforts to fit Sgr A*'s flares with hy- 
drodynamic instabilities find an inclination of roughly 
77° ± 10°, though the systematic uncertainty associated 
with their model is unclear. Fitting a qualitatively dif- 
ferent model to the long-wave l ength observations, a hy- 
drodynamic jet. lMarkoff et a l. (2007) also favor large in- 
clinations (9 > 75°). However, it is not obvious how this 
qualitative agreement should be interpreted, given the 
very different geometries involved. 

Unlike the inclination, our position angle estimate dis- 
agrees significantly with many previous efforts. The pri- 
mary solution is in good agre ement with the 3.5 m m and 
7 mm fits of RIAF models bv lHuang etafl (|2007f ). which 
imply —50° < £ < 30°. While there is considerable 
uncertainty in their estimates (like ours), we note that 
their result is not consistent with our second solution at 
the 1-cr level. In contrast, only our second solution is 
in agreement with analyses of NIR polarization obser- 
vations, assuming a simplisti c orbiting hot-spot model, 
which find 60° < £ < 105° (jMever et al.ll200l . This 
is not particularly surprising given the difficulties in as- 
sessing the relationship between the emitted polarization 
and the underlying plasma geometry. At the 1-cr level, 
none of our solutions are consistent with efforts to fit jet 
models to Sgr A* at 7 mm, which find 80° < £ < 120° 

7 Indeed, as mentioned above, the most likely overall set of pa- 
rameters has 9 = 90° . 
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(|Markoff et al.ll2~007f ). This is also not unexpected given 
the considerably different geometries involved. Between 
all of these (even prior to our estimate), nearly all values 
of £ have been reported, suggesting that significant im- 
provements in theoretical understanding and additional 
observations are required before this question can be set- 
tled. 

Finally, note that our estimat e s imp ly that the X-ray 
feature described in lMuno et al.1 (|2008| ) is not associated 
with a putative jet, but simply another filamentary struc- 
ture in the vicinity of Sgr A*. However, inte restingly, our 
prim ary position angle estimate (and that o flHuang et al.1 
l2007f) is similar to the orientation of the compact, clock- 
wise disk of massiv e stars orbiting the Galactic center 
(|Genzel et alj |2~003). This is somewhat surprising given 
the the large radii (2 x 10 5 GM/c 2 ) of the stellar disk, 
though may be natural if the observed stars are the rem- 
nants of an active period in Sgr A*'s recent past. 

6. CONCLUSIONS 

Despite the sparse u—v coverage, and the existence of 
a detection on only one long baseline, we have been able 
to significantly constrain the possible parameters of an 
accretion flow onto Sgr A*. Within the context of a qual- 
itative RIAF model that fits the observed spectra of Sgr 
A*, we have have been able to substantially constrain the 
orientation of Sgr A*'s spin. The magnitude of this spin 
is less well determined, though the black hole cannot be 
maximally rotating. This result is relatively insensitive 
to the black hole mass. 

We have not ascertained the strength of our con- 
straints' dependence upon the particular RIAF model, 



though our estimates are consistent with earlier efforts 
comparing longer wavelength observations to an alter- 
native RIAF model. However, there are two additional 
reasons to believe that our results will be generic of 
RIAF models. The first is that the underlying physics 
that limits the orientation is the Doppler beaming and 
boosting that is dependent primarily upon the Keplerian 
velocity profile of the accretion flow, a feature that is 
generic among RIAF models. The second is the weak 
dependence upon large variations in mass, implying that 
changing the scale lengths of the accretion model does 
little to alter our results. Thus, we expect the qualita- 
tive form of our constraints to be a generic feature of all 
RIAF models for Sgr A*, and our quantitative results to 
be roughly correct. Of course, should the emission ob- 
served from Sgr A* not arise in an accretion flow, our 
results could be quite different. 

As implied by the right-panel of Fig. [2 additional 
long-baseline observations are sorely needed to unam- 
biguously determine both the applicability of RIAFs to 
Sgr A* an d constrain their parameters. In a compan- 
ion paper, iFish et aTl (|2008t ) we report upon an analy- 
sis of the ability of possible millimeter VLBI arrays to 
do so. Given the number of potential radio telescopes 
that could be added t o the current 1.3 mm VLBI array 
(|Doeleman et al1l2008t ). the prospects for making signif- 
icant progress in model parameter estimation are excel- 
lent. Already it is clear that with the advent of millimeter 
VLBI we are now entering an era of precision black-hole 
accretion physics. 



APPENDIX 
NON-DETECTION PROBABILITY 



Consider an observable y, with an expected value of yo and Gaussian random errors of amplitude Ay. Then, the 
probability density of measuring a value y is 



p(y) 



-(y-y ) 2 /2Ay 2 

V2nAy 



(Al) 



A non-detection, by definition, corresponds to a "measured" value of y below some threshold Y. The probability of 
such an event is simply 



P(<Y) 



p(y)dy = \ + 



y p -(y-y ) 2 /2Ay 2 



2nAy 



-dy = - 

J 2 



1 



erf ( Y -yA 
\V2Ay) 



where 



erf (x) 



■* dt, 



(A2) 



(A3) 



is the standard error-function. As expect ed, when Y = yo, p(< yo) = 1/2 and when Y ^> yo, p(< yo) ~ 1. This 
procedure is identical to that employed by[Kelly (2007), §5.2, and references therein. While in that work the primary 
interest was a careful investigation of the ability of Baycsian techniques in the context of linear regression, we are faced 
with a more general fitting problem. Nevertheless, how non-detections, or "censored data", enter into the likelihood 
is identical, the distinctions arising only later in the analysis. 
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